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ABSTRACT 


In  the  interest  of  developing  more  thorough  methods 
of  analyzing  ship  structures,  an  entire  bottom  section 
extending  from  side  to  side  and  from  one  bulkhead  to 
another  is  analyzed  as  an  orthotropic  shell  of  variable 
curvature.   The  longitudinally  framed  destroyer  type  ship 
is  studied  under  hydrostatic  loading  and  longitudinal 
bending  stress.   Specific  consideration  is  given  to  the 
orthotropy  of  the  structure  and  its  curvature. 

As  a  result  of  the  combined  analysis  of  the  shell 
and  its  associated  stiffener  system,  .partial  differential 
equations  are  developed  which  define  the  behavior  of 
the  stiffened  shell.   The  inherent  complexity  of  the 
developed  equations  precludes  an  analytic  solution.  An 
approximate  solution  could  be  obtained,  however,  by  the 
use  of  difference  equations  and  digital  machine  computation. 

An  application  of  membrane  theory  to  an  equivalent 
orthotropic  shell  provides  a  means  of  determining  direct 
stresses  and  extensional  rigidities.   A  comparative 
study  of  two  existing  destroyer  types  using  developed 
theory  reveals  that  membrane  theory  with  corrective 
iterations  is  of  value  in  preliminary  analysis.   The 
sensitivity  of  the  variation  of  hull  curvature  from  the 
analytic  curvature  is  also  demonstrated.   The  use  of 
the  midship  section  coefficient  as  a  design  parameter 
for  the  structural  designer  should  be  encouraged  and 
exploited.   If  the  midship  section  area  is  kept  constant 
the  curvature  can  be  varied  to  the  advantage  of  the 
structure  with  no  effect  on  speed.   The  relatively  high 
bending  stiffness  in  the  girthwise  direction  of  current 
designs  merits  a  reevaluation  of  transverse  strength 
requirements • 

Structural  design  procedures  employing  orthotropic 
rigidity  coefficients  in  combination  with  either  the 
solution  to  the  above  differential  equations  or  a  modi- 
fied membrane  theory,  have  been  outlined  by  the  author. 
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NOTATION 

b^  -     transverse  web  frame  spacing  (in.) 
b2  s    longitudinal  stringer  spacing  (in.) 

di    =      <L_  =  (     ) » 

d2   =   2-  »  (    > 

^  0 

din  =  21 
3xn 

-   IL* 

p  *  external  load  (lb/ft2) 

q  „  hydrostatic  load  (lb/ft2) 

r  -=■  radius  of  curvature  (ft) 

t  =  shell  thickness  (in.) 

t  ■=  equivalent  shell  thickness  (in.) 

u  -  axial  displacement 

v  -  girthwise  displacement 

w  -  radial  displacement 

x  -  longitudinal  coordinate  (ft.) 

z  -  distance  from  shell  middle  surface 

A.  =  cross  sectional  area  of  a  stiffening  member  (in.  ) 

B  r  beam  of  ship  at  design  waterline  (ft.) 

C     =   dist.  of  centroid  of  stiffening  member  from 
shell  middle  surface  (in.) 

D     -       extensional  rigidity  (lb/in) 


B   =  modulus  of  elasticity  (psi) 

G-   =    shear  modulus  (psi) 

H   »  draft  of  ship  at  design  waterline  (ft.) 

I   -  moment  of  inertia  of  stiffening  member  about 
its  N.  A.  (in**-) 

J   -=  polar  moment  of  inertia  (in*4-) 

K   -  bending  rigidity  coefficient  (lb  -  in) 

L    ■=  length  of  shell  between  bulkheads  (ft.) 

Ls   =  length  of  ship  (ft.) 

M    -  moment  resultant  (lb-ft/ft) 

N    =  force  resultant  (lb/ft) 

Q,         -  transverse  shear  force  resultant  (lb/ft) 

S     =  rigidity  moment  coefficient  (lb-in) 

Vg   __  transverse  vertical  edge  loading  on  shell  (lb/ft) 

VL   -  longitudinal  vertical  edge  loading  on  shell  (lb/ft) 

Y  -  shear  strain 
£   =  linear*  strain 

V  -  Poissons  ratio 

^~   =  stress  intensity  (psi) 

^~M  -  elastic  yield  stress  (psi) 

(j~  -  longitudinal  hull  girder  bending  stress  (psi) 

ff    -  shear  stress  (psi) 

<D    =r  girthwise  or  circumferential  coordinate  (deg) 
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I.   INTRODUCTION 

This  analysis  can  be  considered  as  a  contribution  to 
the  continuing  effort  being  made  to  produce  more  rational 
and  less  empirical  ship  structural  design  procedures.   The 
analysis  of  a  structure  whose  loading  and  response  character- 
istics are  not  readily  defineable  does  indeed  present  a 
challange.   There  exists,  however,  opportunities  for  improve- 
ment in  the  static,  elastic  analysis  of  ship  structures. 
Such  improvements  could  contribute  towards  a  better  determina- 
tion of  adequate  but  not  excessive  scantlings  arranged  in  the 
most  efficient  manner.   Scantlings  are  normally  selected  on  +^e 

basis  of  a  series  of  approximate  analyses.   While  the  effect  of 
continuity  is  often  considered,  the  effects  of  adjoining  members 

and  shell  curvature  are  neglected.   Checks  are  made  on  the 

adequacy  of  stiffeners  and  stiff ener-plating  combinations  on 

a  panel  basis.   This  artificial  method  of  analysis,  while 

adequate  and  sufficiently  approximate  in  the  past,  cannot  be 

cpnsidered  as  satisfactory  for  the  future.   The  knowledge  of 

ship  seaway  loadings  is  increasing.   Advances  are  being  made 

in  dynamic  loading  theory.   If  these  accomplishments  are  eventually 

to  be  used  to  greatest  advantage  in  design,  they  must  be  used 

in  conjunction  with  a  realistic  static  analysis. 

The  structural  analysis  of  a  plate  or  shell  stiffener  system 
can  be  approached  in  general  by  either  of  two  basically  diff- 
erent view  points.   The  first  is  that  of  a  grillage  beam  analysis 
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wherein  the  plating  is  taken  to  act  as  additional  flanges  to 
the  stiffeners  in  an  orthogonal  network.   Essentially  the 
grillage  beam  concept  as  originally  developed  by  Vedeler  vlj* 
relys  on  equating  deflections  at  stiffener  intersections  and 
moment  distribution.   The  solutions  may  be  obtained  by  ortho- 
dox solutions  to  differential  equations,  trigonometric  series, 
energy  methods  or  relaxation  and  numerical  methods  suitable 
fpr  digital  machine  computation.   Variations  in  the  grillage 
beam  technique  and  its  solutions  have  also  been  developed  (2). 
One  need  only  to  consult  contemporary  technical  journals  re- 
lating to  aeronautical  and  civil  engineering  to  find  a  prolific 
assortment  oftechniques  and  methods  for  grillage  beam  solutions 
most  of  which  lead  ultimately  to  machine  computation  (e.g. 
stiffness  matrices). 

The  second  basic  approach  is  to  utilize  plate  or  shell 
tb-eory  for  the  given  structure  in  such  a  manner  that  the 
stiffener  system  effect  is  incorporated  into  the  medium.   This 
ip  accomplished  by  substituting  for  the  actual  shell  with 
stiffeners  (frames  and  longitudinals  J  an  equivalent  shell  with- 
out stiffeners  which  exhibits  different  elastic  properties  in 
the  longitudinal  and  circumferential  directions  (orthotropic 
behavior J .   The  manner  in  which  the  equivalent  orthotropic 
shell  is  determined  will  be  discussed  shortly.   Let  it  suffice 
at  this  point  to  state  that  the  result  of  such  a  substitution 
will  be  to  obtain  a  shell  structure  whose 'general  behavior 

«  Numbers  in  brackets  after  an  author's  name  indicate 
references  listed  in  Appendix  E  of  this  paper. 


under  loading  is  equivalent  to  that  of  the  stiffened  shell. 
The  orthotropic  behavior  ofthe  equivalent  shell  may  be.  called 
constructional  since  it  does  not  result  from  a  property  of  the 
material  but  from  the  manner  in  which  it  is  distributed.   It  is 
interesting  to  note  at  this  point  that  this  orthotropic  approach 
is  diametrically  opposite  to  the  grillage  beam  method  (i.e., the 
grillage  beam  assumes  plating  to  act  as  additional  flange 
material  for  stiffeners  whereas  in  the  orthotropic  shell,  the 
stiff ener  influence  is  distributed  to  the  shell  or  plating) . 
It  is  felt  that  the  orthotropic  shell  concept  has  sufficient 
merit  to  be  the  general' approach  used  in  this  analysis.   Physi- 
cally it  appears  to  be  more  meaningful  in  the  response  and 
behavior  of  the  stiffened  structure.   While  the  task  of  analysis 

of  such  a  structure  is  not  simple,  it  is  at  least  amenable  to  a 
mathematical  investigation.   Orthotropic  shell  theory  consistent 
with  the  structure  being  studied  will  be  utilized. 

Additional  methods  of  shell  or  plate  stiffener  systems 
include  the  concept  that  the  stiffener  system  acts  as  an  elastic 
support  for  the  loading  and  as  such  may  be  treated  by  an  applica- 
tion of  the  beam  on  an  elastic  foundation  method  (3K   Another 
interesting  technique  is  to  disseminate  the  effect  of  one  set 
of  stiffeners  (frames)  into  the  plating  while  other  (longitu- 
dinals) set  is  subjected  to  a  harmonic  analysis  producing  influence 
coefficients  for  deflection  or  bending  moment  for  each 
longitudinal  (I4.)  . 

This  analysis  will  investigate  the  structural  behavior  of 
the  bottom  section  of  a  ship  from  side  to  side  and  from  bulkhead 
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to  bulkhead.   The  ship  type  considered  will  be  a  longitudinally 
framed  naval  warship  in  the  destroyer  catagory.   J&aphasis  will 
be  given  to  shell  curvature,  continuity  of  the  shell  between 
bulkheads,  and  the  stiff ener  characteristics.   Normal  hydro- 
static loading  combined  with  the  overall  longitudinal  bending 
stress  will  be  considered.   Shade  (5)  developed  a  bendine 
theory  for  ships  bottom  structure  utilizing  orthotropic  plate 
theory.   However,  a  flat  bottom  was  assumed  and  no  considera- 
tion was  given  to  hull  curvature. 
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II.  ANALYSIS 


Preliminary 

The  portion  of  the  ships  structure  to  be  analyzed  is 
sjiown  in  Figure  2  .     u &  m 


r 

**9 

1       1        ,      Tl    - 

I1       1       1       1        1       l| 

'i                    i        i       n 

_ 

/  Keel 


Fig.  2  Identification  of  Section  Analyzed 

The  boundary  conditions  for  this  section  are  discussed 
in  Appendix  C  .   In  brief  they  are  approximated  as  follows : 

The  degree  of  fixity  can  be  assumed  to  be  100%  clamped  at  the 
junction  of  the  shell  and  bulkhead  and  simply  supported  at  the 
junction  of  the  curved  shell  to  the  vertical  or  flat  portion 

of  the  ships  side.   While  these  assumed  conditions  are  not  in 

reality  present,  a  determination  of  the  elastic  degree  of 

fjLxity  is  outside  the  realm  of  this  analysis.   The  postulated 

conditions  will  be  sufficiently  accurate  for  approximation 

purposes. 

The  inclusion  of  the  huj.1  curvature  in  the  analysis  can  be 
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effected  by  the  selection  of  an  analytical  expression  for  the 

radius  of  curvation  (r),  as  a  function  of  girth  angle  (0) . 

The  selection  of  the  expression  for  the  radius  of  curvature 
must  be  uncomplicated  enough  so  that  it  will  not  impede  the 
mathematical  operations  involved  in  the  analysis.   At  the 

same  time  it  must  sufficiently  describe  the  hull  form  for 
any  given  B/H  ratio.   A  consideration  of  pertinent  hull  forms 
and  analytical  expressions  is  made  in  Appendix  A  .   The  result- 
ing analytical  expression  satisfying  the  above  requirements 
has  been  found  to  be  r  -  2Hcos0;  the  radius  of  curvature  which 
describes  a  cycloid  of  depth  H  and  span  if  H.   One  of  the  sig- 
nificant contentions  of  this  thesis  is  that  a  ship  type  such 
as  a  destroyer  has  a  hull  form  which  cannot  be  treated  with 
f^Lat  bottom  or  flat  plate  assumptions.   Of  course  if  a  restric- 
tion is  made  to  only  a  small  panel  of  the  structure,  a  flat 
surface  is  more  closely  approached.   This  type  of  panel  analysis, 
however,  is  a  rather  limited  view  of  the  structure.©  behavior 
and  disallows  for  the  continuity  of  deformation.   Moreover,  not 
only  is  significant  curvature  present  but  it  is  variable  around 

the  circumference  of  the  hull.   An  assumption  of  constant  cur- 
vature or  a  circular  cylindrical  shell  would  greatly  facilitate 
the  mathematical  analysis  but  would  not  yield  results  which 
are  representative  of  the  structure.   It  is  felt  that  the 
variable  curvature  selected  represents  sufficiently  the  type 
of  hull  forms  being  considered  (see  Figs . 19 "23) . 
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Isotropic  Shell  Analysis 

To  begin  the  analysis  the  elastic  laws  for  the  variable 

curvature  isotropic  shell  (ships  hull  without  stiff ener  system) 

will  be  developed.   The  shell  element  and  associated  stress 

system  are  indicated  in  Pig.  J>, 

Equilibrium  of  forces  provide  the  following  relations: 

x  dir: 

"£>   (Nx)    dxrd0      Q    (N^x)    d0dx  f  pxrd0dx  ^  0 


c>  x  ^TF 


or 


(1) 


+     Px    =    ° 


Snx  i^x 

3x~    "—  dT 

0  dir: 

cMN0)    d0dx    3    (NY^)    dxrd^T  -    Q0dxdJ2[     -      (40dxd0 

"  §  LSri)    d0dxd0  p^rd^dx       -    0 

'   (2) 


or 


i  3ty      ^  -  i  <fy  ±  P0   ^  o 

p    Q^  6x  r 


-9. 


r  dir: 


-  2N0dxdGi        -      d  Uy)dxrd0     -     3   (<^)d0dx  +  p~rdxrd0  =  0 

2  ^1T  a? 


or 


i    N0  +   9_^k    +  i    ^ 


Pr     -     0 


3x  r    c>0 

Equilibrium  of  moments   give     — 
about  p^   axis : 


-  c^K^)    dxrd0      -     C>  (M^)d0dx  +  Q^dxrdjfl    *    0 
c>x  3  ^ 


or 


0- 

3   (Mx)dxrd0  ,     3  (Mpx)d0dx     -   Qxdxrd0  =  0 
5x  dF~ 


or 


9   Mx        1  3_M0X  -   Qx  =  0 
"£x  r  5^ 


•  10- 


(3) 


c>x  r    o)0  (I|-) 


about  pM   axis : 


(5) 


About  p      axis : 


-Nxo(rd0dx  4    N0xdxrd0  -   2   %xdxd0  -   0 

2 


or 


(6) 


N..^  -   N 


<jt  "   >c+  i  n0x 


1  MM„  =   0 

r 


Elimination  of  Q^  and  Q^  by  use  of  eq.  (I4.)  and  -(5)  gives  the 
'following  equations  of  equilibrium. 

3-Nx+  i^V     +    Px    -    °  (7a) 

Ox  r  dp 


1  djfy    l_\y  -  1  ^  -  1  ^y  f  P/   =  0 


r  d^         dx 


•2  c)   (/ 


(7b) 


1  %^  ,  1  ^i.M0x4  i 


^ 


xL 


£>xd0 


30Jk       dj# 


k 


•p„  =  0         17c) 


r 


(7d) 


At  this  point  there  are  k   equations  and  8  unknowns.   A  study 
of  the  shell  deformation  will  provide  additional  relationships. 
In  the  analysis  of  deformations  the  following  assumptions 

are  made : 

(1)   All  points  lying  on  one  normal  to  the  middle  sur- 
face before  deformation  do  the  same  after  deformation  (i.e.,  neg- 
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lect  deformation  due  to  Q^.,  Qu   transverse  shears), 

(2)  For  all  kinematic  relations,  the  dist.  Z   of  a 
ppint  from  the  middle  surface  may  be  considered  as  unaffected 
by  the  deformation  of  the shell,  but  that  for  all  considerations 
of  the  stress  system,  the  stress  0" z  in  the  2  direction  may  be 
considered  negligible  compared  to  Q~"x  and  $~d. 

(3)  All  displacements  are  small,  i.e.,  negligible 
compared  with  radii  of  curvature  of  the  middle  surface  and 
their  first  derivatives  are  negligible  compared  to  unity 
(linear  assumption). 


Seobun  cp  =  co<j=,i. 


jec-Uolr\,      X  =  Cgn%\: 


Fig,  k  Displacement  of  Point  A 

U,  V,  W  refer  to  displacements  in  middle  surface  of  shell 

in  directions  x,  0,  and  z  respectively. 

UA>  VA»  WA  refer  to  displacements  of  an  arbitrary  point  A  at 
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a  point  Z  from  the  middle  surface 
from  Pig.  l\.    : 

UA=  U  -  Z  d   W 


vA  -.  V  (^z) 


z\  -  Z  d  W 


tfj^  -  W     (assumption  2) 

The  strains  at  point  A  can  be  expressed  by  the  following 
relations  from  Flugge  (6). 

_^*A 

6   -  L__4lA    -  WaN 

0   r  +  zU 


r« 


0 

^)VA    1      >A 


x0   Sx     r+z  c)^ 


The  strains  may  now  be  expressed  in  terms  of  middle 
surface  displacement  by  using  point  A  displacement  relations 
^u   -z  c)2w 
x 


£_- 


/-  1 1        7)        I  v(lt£    )     -     Z^W    |   |WL 

y     r+z  \  "bp  L      r        r  STPTJ      J 

V  ^_       v(l+   zl   -   z     £wl      1 r  bjx  -   z    c)   w 

x-p      ^x  L  r/        r     ^J       r+z   L  d  0  B0d* 
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A  simplification  can  be  made  at  this  point  by  considering 

the  shell  thickness  to  be  very  small  in  relation  to  the.  radius 

at  any  point.   In  otherwords  z  is  negligible  compared  to  1. 

r 
or 

(1  +  z/r)  a.  1 

_1 _  1 

r  t  z   ""*  r 


The  strain  relations  now  become 


2  u  -  z  2>  w 


3jC 


6 


H 


0 


4 


1 
r 


_^0 


-  z_  o)  w\  i  w 
r  3>0  J 


-  z 


^w 


Carrying  out  the  indicated  operations  and  remembering  that 
r  is  a  function  of  ft   (r  -  2HcosjZf) ,  we  obtain  the  following 
relations  when  the  notations   B__  _(   )'  and  ^_   (   )  are 

dx  ~~      ■        aw "' 

employed. 


£       -    u! 

X 


z  w' 


(8a) 


C   ,  _  v  -   z   (w  tan  0  f  w)   w 
J*  "  ?     r2  r 


(8b) 


r. 


0 


u 
r 


z_ 
r 


w 


(8c) 
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The  stress  resultants  of  Fig.  3>  are  by  definition 
t/2 

Nx  =  J     <TX  (1  *  z/r)  dz    »  N^  J  <T^  dz 

-t/2  -t/2 

t/2  t/2 

Nxtf  =  J^Vi  +  z/r)  dz    ,  N0*  -  J^  ^A  dz 

-t/2  -t/2 

(9) 

t/2  t/2 

^  -  m       C    (Tx(l*h  z/r)  zdz    ,  M^r-   j   ^zdz 

-t/2  -t/2 


t/2  t/2 

M  ^  .-  -     T/rx0(1  +  Z/P)  ZdZ    '        ^x"'   yr^xzdz 
-t/2  -t/2 


The  terms  (1  +■  z/r)  goes  to  1  however  by  the  previous 
assumption  of  negligible  z/r  ratio. 

These  stress  resultants  can  now  be  found  with  the  aid 
of  the  strain  relations  (  2    ) • 


■45- 


since 


(Tx-r^~2(£x^  ty 


(\o) 

1-1) 


1x0   *  -S 0  xCf 

2(l+-»)     P 


After  substitution  of  {  \o  )    and  (  8  )  in  (  9  )  tne 
following  stress  resultant  relations  are  obtained  where 


c  =  liji      is  the  extensional  rigidity  and 


K-  _J*t _   is  the  flexural  rigidity  or  bending 

12(l-^>2) 


stiffness. 


N  D      (ru'tVv      +   V  w) 

x    =    — 
r 

N<*   -     £     (v  +  w    -f  -v>  ru«)  (H>) 


Nx0    '      N)ZJx    -    D    (1  "V)         (^  *"  rvl) 

2r 
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M^.  _  K   (r2  w"  4--VW  \  -O   tan  <f)   w) 


M. 


^ 


K   (w  +  tan  0  w  +"9  r^  w") 


(11  con't) 


Mx0  -  V*  ~-  K  (— >   W 


i* 


These  relations  represent  the  stress  resultants  of 

an  Isotropic  cylindrical  shell  with  variable  curvature 

r  -  2HcosJ^  .   (Note  that  the  shell  thickness  has  not  been 

fixed  yet  but  can  still  be  a  function  of  0).   These  relations 
exist  in  the  unstiffened  ships  hull. 

We  consider  now  the  gridwork  of  frames  and  longitudinal 

stringers  only^ in  Figure  5 


k>, 


sS 


SUcA.io'J   tf,  C*M^"t 


S^oA»0*>   X"  =  Oe>~5"t 


Fig.  5  Element  of  a  Gridwork  Shell 
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N2  =   KA.0  £  ,      where  L  <i     -  strain  at  centroid 


then  N, 


0   =  ^2      ^  ^  C( 


Using  the  strain  relations  of  (eq.8)with  z  =  CW 

C0  |ii 


H0  =  ^t 


V  -  C^  (w  -+-  tan  0  w)  w 


r   r' 


similarly 


N 


EAX6 


x^  x 


EAx  (u»  -  CxwH) 


^    M2  "  N2Ctf 


EI0"  -   EA0C0  £  0 
bXr2       bn 


where  w  _  change  in  curvature 
~2  " 


M. 


0  -   ^0 
b-^r' 


2,*< 


EIaw   -   EA^C^  (v/r+  w/r)   EA^C^  (w  -htan^w) 


4- 


b^r' 


similarly  for  Mx, 


2„n 


^  -  Ml  "  Nlcx    EIxw"  "  EAxCxU1     EAxCx^w 


bp       bp 
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Mx0  = 


M 


12 


^x  = 


M21 


twist  of  a  stringer  =.  6]_  ^  M12      ^±   w"  ^ 


GJ- 


6lbl  _  M12bl        Mx0b1b2  bxw 


N    = 


GJV  GrJ. 


G<3X     w 


b2       r 


similarly,    M0X    GJ^ 


w 


npte  that  M0X    J0   b2 
*W    Jx   bl 

Summarizing  the  gridwork  stress  resultants  of  a  variable 
curvature  cylindrical  shell  form  (r  =  2Hcos^j .   — 


Nx    EAX  (u»  -  Cxw") 


*ty    EA^ 


-  C0  (w  4-  tan)6w)  4.  w/r 


~2~ 
r 


Mj  _  Elxw*   -   EAxCxU'     EAXC 


2  w« 


x 


-e 


(12) 
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Mp  ^     El^w      -      EA^C0    (v/r  f  w/r)        EA^2    (w  +  tanjtfw) 

—      4_  __ 


Mxtf    = 

2fx 

»•'■ 

b2 

r 

M0X    _ 

GJ0 

./ 

w 

bl 

r 

(12  con't) 


NOTE:   Nx0  and  N^x  could  also  be  determined  but  it  will  be 
seen  shortly  that  they  will  not  be  necessary. 

These  relations  exist  in  the  ships  hull  stiffners. 
Superposing  the  stress  resultants  of  the  isotropic  shell 
with  those  of  the  gridwork  gives  final  stress  resultants  of 
the  stiffened  variable  curvature  shell  (or  the  stiffened  ships 
hull ) . 

N   -  Dxu«  ,   D^>  (v  +  w)  -  SxW" 
4-  


]<p       D^(v  +  w)  -4-  D^u»  -  S^  (w  4-  tan0w) 


r 


(13) 


Mx  =  Kxw"  .   Ko  (w  +  tan^w)  -  Sxu» 
3~ 


M^  _  K0  w   +  k   |w"+  tan0w)  -  S^  (v  f  w 
r2     V  *    T^'        ~ 
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N 


x0  _   Dxjtf  («  +  rv») 


N0x  .   Dx0  (*  4-  rv»)  ^  K0X  w 


. /  (13  con't) 


r2 


*fyx    ^x  »' ' 


Mxtf  _  *x0  w#  ' 


Where  the  orthotropic  rigidity  coeffecients  are  defined  as 


Dx 

Et     EAY 
=  2+— * 

D 

Etv 

l-^2 

Kx 

Et3 

4- 


12(1  -^2)  b2 


D^  _  Et     EAfl 
l--p2  *  bx 

2(1+ V  ) 

K0    Et3        E(l0 

=  ,  ■■         i 

A0c/) 

12(1  -v*) 

bl 

(14) 

*  2. 

E! 


K      ^  Et3v  ,  Sx     _      EAXCX  ,         S^     ^      EA0C0 

^  *    12(1--/)  "bi"  bl 


Kx0  _      Et3  GJX  Kjfe  Et3  GJp 

12(HV   )  b2  12(1 +  V>    )  b-L 
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Note  that  the  Nxj  and  N0X  represent  shear  carried  by 
the  shell  alone.   The  gridwork  can  only  carry  this  shear 
by  a  bending  of  its  stiffeners,  a  bending  which  becomes 
impossible  when  the  shell  is  joined  to  the  stiffeners. 

Introducing,  now,  these  new  stress  resultants  for 
the  stiffened  ships  hull,  into  the  equations  of  equilibrium 
(eqs.  7),  one  obtains  the  following  differential  equations. 


.u"   D^  (v»  4.  w«)  -  Sxw«  "    1 


r  90 


D. 


:(j(u  +Vj 


iiL 

r    clifi 


D^  (♦+„)  +  Dvu- 


S^(w  4  tan0w) 


\ 


Dx0(u«  4  rv»)  -  K^  w» 


Sj  ( v  4-  w ) 


4  P0 


PrA  =  0 


D^(v4-w)   D^  u'  -  £U(w  +  tan0w)  +  Kxw«  "  »+  Kv  (w"  +  tan0w" ) 


-  Sxu»'»    1  S  (K0X 


r  d0\  r 


(¥ "7 


+ 


r  ^   r       rc)^Lr2  V  r2' 


_  °0 


S«<  (v  +  w) 


p    -        0 
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After  carrying  out  the  indicated  operations,  rearranging 
terms  and  using  partial  operator  notation  as  follows, 


and 


din  =  si 

a  x 


d2  -_  9_     .  (      ) 

d2n_    ^n 
^0n 


the  differential   equations   assume   the   form, 

rir3dl  +  G2rd22    "*"   G2rtan0d2~|    u    4      |  G3r2dld2]     v 

Fl4-p2dl  "   ^r3dl3   *"     °6dld22  v   2G6tanQf  did2~lw 


(15a) 


■r^p 


x 


|C3r3d1d2       u    +   |G7r2d22  4    c8^d22  4    C2r^dl2  +"   V^   tan(Z*  d2 
+    CQr  tan0  dgl     v    4-    r-c9d32   "   c8rd23   "   c10r2dld2 +  c7r2d; 
-3CQr   tan20  d2  -   C13d2  -  3C1;Ltan20  d2   -   Cgr  tan0  d£2 


(15b) 


-2C9  tan0  d22  -  C11  tanp  d22  -f  C7  r  tan0  +  Cgr  tan0  w-  -rT?0 
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r-C^d134    C^r3d1~1       u   4-  P-C8rd23   -   3C8r  tan0  <1* 
■\  C7r2d2   -   3C8r   tan2fi  d2  -   Cgrd^l      v  4. 

rc12r^d1Il-  ^   C13r2d12d22  +   C^  4    C13r2tan0  d^ 
+   Cx^tan0  d23   -   2C8rd22  +   8C9tan20  d22  +    2C  d£2 


(15c) 


4  2Cl:Ld22  4  7Cl;L  tan20  d22  -  ^Cgrtan0  d2  4-  9Cxltan0  d2 
4-  l^Cn  tan30  d2  4-  C7r2  -  3Cgr  tan20  -  CgrH  w  =  r^pr 


The  constant  part  of  the  coefficients  are  identified  as 


C1=DX 


c8  =  S 


0 


C2  -  Dx0 


C-,  = 


D-*+  Dx0 


\=^ 


c5  -sx 


C6  = 


C7  = 


K 


0x 


CQ  - 


9  =  ^ 
;10  =  K%>  +  Kx.0 


C12  -  Kx 

01U  =  %  4  5i 


Equations  (15,  a,  b,  c)  comprise  a  set  of  3  partial  differential 
equations  in  3  unknowns,  (u,v,w).   Their  solution  as  a  function 
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of  radius  and  rigidity  coefficients  will  completely  define 
the  structural  behavior  of  a  stiffened  cylindrical  shell 

(of  cycloid  cross  section)  for  given  loading  and  boundary 
conditions.   The  rigidity  coefficients  were  made  constant 
as  a  necessary  simplification.   It  will  be  seen  later,  how- 
ever, that  many  are  in  actual  practice  essentially  constant 
over  considerable  portions  of  the  shell  section.   Additional 
comments  concerning  these  differential  equations  will  be 
delayed  until  the  Discussion  of  Results  section. 
Membrane  Analysis 

Since  the  solution -of  the  developed  differential 

equations  is  not  readily  available,  one  must  approach  the 
shell  problem  with  other  and  more  immediately  useful  tools 
of  analysis.   It  must  be  stated  at  this  point,  however,  that 
the  analytical  approach  by  means  of  the  differential  equations 

of  the  problem  is  the  only  means  by  which  a  complete  and 
relatively  exact  solution  may  be  obtained.   It  is  unfortunate 

that  their  complex  nature  is  currently  limiting  their  use 
by  shell  designers.   The  alternate  methods  employed  to 

obtain  solutions  include  the  beam  method,  iterative  calcula- 
tion techniques,  membrane  theory,  and  other  unique  approxima- 
tion approaches.   Lundgren  (7)  describes  rather  completely 
the  various  methods  of  shell  analysis  currently  in  use. 

The  applicability  or  the  relative  desireability  of  the 
alternate  methods  is  determined  basically  by  the  shell  dimen- 
sions.  The  ratio  of  the  shell  length  to  the  shell  span  (or 
beam),  the  L/B  ratio,  normally  will  govern  the  technique.   The 
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application  of  membrane  theory  (i.e.,  the  consideration  of 
direct  stresses  only;  the  support  of  the  load  by  tension, 
cpmpression  or  shear  forces  only)  is  considered  good  for 
short  shells  or  shells  whose  L/B  4C.  1.   Lundgren  (7)  states 
that  the  membrane  condition  of  stress  is  a  very  good 
approximation  for  the  short  shell  since  the  elastic  support 
between  end  sections  (or  bulkheads)  is  quite  stiff  and 
loading  is  balanced  by  N    N,    and  N  .  stresses.   The 
ASCE  Manual  of  Engineering  Practice  No.  31  recommends  a 
superposition  of  membrane  stresses  and  those  created  by 
boundary  conditions.   However,  on  the  matter  of  boundary 

conditions,  the  manual  also  states  that  for  a  L/B  <  .83 
the  transverse  (Bulkhead)  end  effects  damp  rapidly  and  the 
effects  of  the  longitudinal  edges  are  felt  only  in  a  small 
area. 

The  previous  discussion  indicates  that  a  membrane 
analysis  of  a  ships  bottom  section  between  bulkheads  is  a 
reasonable  approach  when  one  realizes  that  the  L/B  range 
for  the  ships  or  ship  type  under  consideration  is  as  indica- 
ted in  Table  1. 


L/B 


DE  1033 

DDG  952 

DLG  9 

DLG(N) 

CG(N) 

CA  13U- 

0.775 

0.85 

0.81 

0.965 

0.56 

0.56 

L  -  length  between  bulkheads  (  \n  $  Ye<\\oS) 
B  =  beam  at  waterline 


TABLE  1   Length  to  Span  Hatio  of  Ships  Considered 
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The  membrane  stresses  (lb/ft)  in  the  shell  element  can 
be  obtained  by  setting  the  shell  moments  equal  to  zero  in 
the  equations  of  equilibrium  (eqse  7). 

d  x         r   30" 

^>Nx  =  -Px  "  I     ^_V 
5x  r    c)    tf 

integrating  the   latter  two   of   these, 

Nx0  =    "       f(P0+    X    &*#)    dx   +fl(^ 


N 
x 


=    -    ]    (Px+   1^M\     ta   -^f2(^) 


r      c>  0 


for  px  —  0  and  p^  and  pr  independent  of  x.   (These  are  correct 
conditions  for  ships  hull  loading.   The  longitudinal  stress 
resulting  from  the  hull  girder  action  and  determined  from 
the  longitudinal  strength  calculation  will  be  accounted  for 
shortly  but  will  not  be  considered  as  a  px  loading  at  this 
point) . 

Nx0  =  "  (P0  +  1  ,d  Nfl\  x  4-  fi(0)  "=  «xF(0)  +    f£(0) 


U0\ 


Nx  =   x2  dF(ffl)  -  x  dfx(0)    f2(0) 
2  r   d(j)  r   d0 
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F(0)  is  determined  by  loading  and  shell  shape 
(curvature)  but  f;j_(0)  and  f2(0)  will  be  determined  by 
bpundary  conditions  at  the  shell  ends  (bulkheads). 

Consider  first  the  case  of  simple  support  end 
(bulkhead)  conditions. 


At  X  -  t   L/2  , 


Nx  =  0 


then, 


Fig.  6   Shell  on  Simple  Supports 


L2  dF(0)  - 


O- 


8r   d0 


L   dfx(0)    f2(0) 
2r  '   d0 


L2   dF(0)     L   dfL.(.0)    f2(0) 

O^  +  - + 

8r     64  2r       d0 


which  gives, 

fx(0)  -  o        f2(0)-    -  l2   mp) 

8r   d0 
and  therefore, 


Nx0  =   -  *FW 


Nx  =  -  1  {IT  -   kx?)      dF(0) 


8r 


d.g* 


■28- 


Now  consider  shell  end  or  bulkhead  conditions  for  the 
shell  continuous  over  the  supports  and  assume  100$  clamping. 
As  an  aid  in  determining  the  stress  Nx  over  the  supports, 
one  can  make  use  of  the  relationship  given  in  ASCE  Manual 
No.  31:   the  ratio  of  the  stress  (in  a  shell  over  continuous 
supports)  at  the  point  of  support  to  the  stress  in  simply 
supported  shell  at  the  midspan  is  equal  to  ratio  of  the 
moment  in  the  continuous  shell  at  the  support  to  the  moment 
in  the  simply  supported  shell  at  the  midspan. 
OV- 


^E 

To 


M 


E 


M, 


and,  for  a  strip  of  shell, 
Mg  =   -1/12  wL2 


0 


xE   _ 


-  2/3 


M0  = 


1/8  wl/ 


(- 


xO 


-y-^i 


NxE  ~  -  2/3  N 


xO 


and  N  _  was  found  previously  to  be  -  L2   dF(0) 

8r  "d^ 


Fig.  7   Comparison  of  Continuous  & 
Simple  Supports  for  Shell 


N  at  support  U  -iL/2)  is  +  L2   dF(0) 


continuous  over  supports. 


At  X=      £   L/2 


Nx  _       j£        dF()6) 


12r 


d0 


12r  &fti 


for   a   shell 


Hi.  ~  I  -  L>i 


y 


fro 


-w^ 


frfi/i 


Fig.   8      Shell   Continuous  over  Supports 
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or 


L2   dP(0)     L2   dF(0f)   -   L   df^p)        f2(0) 
12r  d0      8r    dgf       2r    d0 


L2   dF((jO     L2  dF(0)    L    dfi(0)   f£(0) 

12r    d0      8r  d0      2r     d0 

which  gives, 

?!  (0)  -  0 

f2  (0)  -  "  -i£  gP(j^J 

21^  d0 


then, 


Nx0     -     -   xF^) 

Nx      =      -    (L2   -    12x2) 
214.  r 

dF(0) 
d0 

N0      =      Pr  r 

(16) 


where, 

F(0)=  (P0+  1  d^) 
r  d  0 

Equation  16  represent  the  general  form  of  the  membrane 
stresses  in  the  ships  bottom  section  of  radus  of  curvature 
r. 

These  stresses  will  now  be  developed  in  a  more  final 
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form  by  considering  the  hydrostatic  loading  and  the 
cycloidal  cross  section  cylinder. 


Fig.  9  Loading  on  Ships  Hull 

o^  ■=  Y  (^o  +  ^) 

Pr»  *   -(<1  ~  W  COS  0) 


Let  q0  =  V  hQ  where  hQ  -  an 
arbitrary  design  head  to  account 
for  wave  forces  on  freeboard  surface 
of  hull, 
z  =  H/2  (1+cos20)  from  appendix  (^  A  ) 

w  =   internal  loading  of  shell 
due  to  machinery  weights  (lb/ft  ) . 


Pr  =  -V£to<)  +  H/2(l  +  cos20)]  +  w  cos  $ 
P0  =  -  w  sin  0 

?**      ° 

N0    "    )    "V    Th0  *"  H/2(l+cos20)]    4-    w   cos  pi     2H  cos  j0 

N^     -         -  V  [    2H  hQ   cos  0  -f     2H2   cos3  f\  4  2wH  cos20 

and  since   F(0)  -  p^  ,  1   d  N0 

*   d  < 


(Ha.) 


F(0)  =  y  \hQ   tan0  4.  3H  sin20J   -  3w  sin0 
Nx  and  Nx(*  are  then  found  by  applying  equation  (16). 

Nx0  =  "  X  T^  (h°  tan0  +  3/2H  sin20>  "  3w  sin0J 
N    L2-  12x2  T3w  -  )f  (hQ  sec30  4.  3H_cos20)  1 

x        i^h-  L  cos^    'J 


(17b) 
(17c) 
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Equation  (17)  are  the  membrane  stresses  for  the  ships 
bottom  section.   With  the  aid  of  these  relations,  the  necessary- 
shell  thickness  (that  is,  the  equivalent  shell  thickness  for  a 
hull  with  no  stiffeners)  can  be  found  as  it  varies  around  the 
girth.   From  this  equivalent  thickness,  extensional  rigidity 

coefficients  (Dx  ,  D^  ,  Da   ,  D^  )  may  be  computed.   For  the 

E~   E~~    E    E~~ 
single  equivalent  thickness  derived  there  results  the  fact 

that  Dx  Da   since  they  both  equal  te  or  l«lte  see  appendix  D. 

i   e  i  — %>2 

However,  in  practice  these  coefficients  are  very  nearly  equal 
as  can  be  seen  from  Figure  26. 

12  te 


si=  a 


12  te  (18) 


ru  _  *** 


4-- 


12  te 


Figure  10  Direct  Stresses  on  Equivalent: 
Shell  Element 


Where  5]_  !a"  tne  null  girder  longitudinal  bending  stress  as 
determined  from  the  longitudinal  strength  calculation. 

Using  the  maximum  shear  failure  theory  and  in  particular 
the  Tresca  version,  


(T^  JVmajL       l/(Tx  "(fa    \±    T*P  or 


<Ty2  -Cx2  -  2(Txr^^>24  ^rfp 


(19) 
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Substituting  equation  (18)  in  equation  (19)  and 

solving  for  te  gives 

,2 


t   ^  -Bf.1  \)B  ■  IjAC 


2A  (20) 

where  A  ■«.  ll&   (  (fi.2  -  (Ty2) 

B  -  21,  ^Tl(Nx  -  H^) 

C  =   Nx2  -  2NXN0  +  N2^  -V  ^Nx20 
It  now  becomes  necessary  to  make  a  judgement  as  to  the 


treatment  of  ^~,  and  ^~  . 

(J~  will  be  assumed 

1.65  will  reduce  its  value  to  20,000  psi. 

.11  be  approximated  by  Tobins  r 

y/L  where  L<,  =  length  of  ship  in  ft 


(J~  will  be  assumed  as  33,000  psi  and  a  safety  factor  of 


(J~  will  be  approximated  by  Tobins  relation  that 


M 1  max  =  ^  s 

(p,  r  Ton/in2 

If  it  is  desired  to  have  ^H"  as  a  function  of  ship  depth 

(or  girth),  with  <H  -  Q~_     at  the  keel  (6  -   0)  and 

1     1  max 

(u  _Xl  max   at  the  waterline  or  N.A.   (0  -  ^2)  (to  allow 

for  the  ship  being  in  the  heeled  condition), 

then  <y  _  y*/    Ls      (5  4-3  cos20)  (21) 

"8T— 
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The  membrane  analysis  has  now  been  sufficiently  developed 
to  make  a  test  of  its  applicability  to  specific  examples.   For 
a  given  ship  (or  design),  the  equivalent  shell  thickness 
variation  can  be  found  by  means  of  equation  (20)  which  requires 
a  determination  of  membrane  stresses  equation  (17)  and  hull 
girder  bending  stress  estimate  equation  (21).   The  extensional 
rigidity  coefficients  are  then  readily  calculated  since  the 


equivalent  shell  has  no  stiffeners  and  D  _  D 


£,J*   */./£ 


E    E   1- 1)"2 

The  calculated  values  may  then  be  compared  to  the  actual 
values  on  the  ship  selected  for  study.   It  should  be  noted 
that  in  the  calculation  of  the  membrane  stresses,  the  longi- 
tudinal section  (X  —  constant)  selected  for  stress  determina- 
tion was  the  section  at  the  bulkhead.   This  was  a  logical 
design  decision  since  a  close  examination  of  equation  (20) 
(or  a  few  comparative  calculations)  will  indicate  that  analysis 
of  this  section  results  in  the  greatest  te  for  any  given  (f 
value  along  the  shell  length.   In  addition,  the  value  of  w  or 
internal  loading  was  assumed  equal  to  zero.   Not  only  will 
this  be  the  more  severe  loading  condition,  but  it  is  also  the 
one  which  current  design  (series  of  flat  panel  analyses) 
employs.   This  assumption  will  assist  in  making  the  comparison 
more  valid.   The  retention  of  this  quantity  in  the  general 
theory,  however,  is  recommended.   Additional  comments  on 
this  subject  will  be  made  in  the  discussion  section. 
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III.   RESULTS 

The  analysis  of  the  orthotropic  shell  by  superposition 
of  the  isotropic  shell  and  a  stiffener  gridwork  (with  all 
effects  of  bending  and  torsional  moments  in  the  shell  con- 
sidered) resulted  in  a  set  of  3  partial  differential 
equations  in  3  unknowns  (u,  v,  w) .   They  are  equations  (15) 
on  page  2Z  .   l>ue  to  their  lengthy  nature  they  will  not  be 
repeated  here.   The  coefficients  in  the  partial  differential 
equations  are  variable  in  that  they  contain  the  radius  of 
curvature  in  powers  up  to  the  fourth.   Equation  (15c)  is  of 
the  fourth  order  while  15a  &  b  are  third  order.   Equation  (15a) 
could  be  made  homogeneous  by  setting  px  -  0  but  equations 
1,5b,  c  must  remain  non- homogeneous  •   All  three  partial 
differential  equations  are  linear. 

The  membrane  analysis  produced  more  tangible  results 
in  the  form  of  expressions  for  the  membrane  stresses  as  a 
function  of:   girth  angle  (jJ),  position  along  shell  length  (X), 
and  ship  characteristics  such  as  ship  length  (Ls),  draft  (H), 
compartment  internal  loading  (w),  distance  between  bulkheads 
(L),  and  a  design  head  on  freeboard  plating  (ho). 

Membrane  stresses  in  ships  bottom  section. # 

Nx  _  L2  -  12x2  p3w  -  Who  sec3jZ)  -f  3H  cos2^"7 
[j^H     L"  v  cqs0'J 

N^  -  -  V  (2H  h0  cos0  ■+  2H2cos30)  •+  2w  Hcos2^ 

Nxd  =   "  x  CY  (hotan£>  f  3/2  H  sln2fi)    -  3w  slnfj 
*  A  bottom  section  constructed  with  no  stiffeners  but  an 
equivalent  thickness  (te)  to  produce  the  same  extensional 
rigidity  as  one  with  stiffeners. 
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Figure II  illustrates  the  membrane  stress  distribution 
along  the  length  of  the  shell  (a  section  $  -  const). 

Figure  \X   represents  the  circumferential  distribution 
of  membrane  stresses  in  JJE  1033  and  DD   931  hulls. (a  Section  x  const) 

The  derived  expression  for  the  shell  equivalent  thick- 
ness variation  utilizes  the  membrane  stresses  along  with 
the  hull  girder  bending  stress  (  (T±)    and  a  yield  or 
allowable  design  stress  (  <Ty) • 


+- 


-B-  \|B^  -  L|AC 


inches 


2  A 


where 


A  ^= 
B  - 
C  = 


iW  (  (Ti2  -   (Ty2) 

2k     (^  (Nx  -  Njj) 


<T,  > 


<Tj   ,    -   Psi 


membrane  stresses  -  lb/ft 


Nx2  "  2NxN0  +  ^0  *  ^4 


The  extensional  rigidity  coefficients  are  then  computed 

l.lte        for  -v?=  0.3 


by  Dx  ^  ^  _  te 


E    E    1-  v> 


D 


4- 


2(1+V) 


V  - 


'e 


.38t, 


.33t 


e 


1-V 


Using  the  DE   1033  and  DD  931  hulls,  these  coefficients 
have  been  calculated  from  the  developed  expressions  above  and 
compared  to  those  calculated  from  the  actual  hull  structure. 
The  results  are  plotted  in  Figure  s  11).,  15, 16, 17  e 
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Figure  11  Longitudinal  Distribution  of 
Direct  Stresses  in  Ships  Hull 
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Figure  11+  Extensional  Rigidity  Comparison 

for  DE  1033 
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Figure   If?     Extensional  Rigidity  Comparison 
V  \  for  DE  1033 
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Figure  16  Extensional  Rigidity  Comparison 
for  DP  931 


£■ 


CoOSfTud^fc. 


^  j  ^      CALCULATED 


— f — 


3o 


Ao 


So 


(P  ^± 


9&| 


#' 


Go 


1° 


So 


Figure   IT     Extensional  Rigidity  Comparison 

for  DP  931 


t) 


I2? 


-CMCuLKTE^ 


<P 


•kl< 


IV.   JjlSCUSSlUi^  OF  RESULTS 

The  inherent  complexity  of  the  developed  set  of 
partial  differential  equations  precludes  a  simple  solution. 
An  analytical  solution  appears  to  be  a  formidable  challange. 
Perhaps  the  most  promising  course  of  action  would  be  to 
apply  the  technique  of  difference  equations  and  with  the 
aid  of  a  digital  computer  obtain  an  approximate  solution. 
The  complex  nature  of  the  equations  should  not  discourage 
efforts  to  solve  them,  for  in  their  solution  lies  the 
complete  behavior  of  the  stiffened  shell.   This  last  state- 
ment needs  qualification  in  that  the  behavior  does  not 
include  the  effects  of  concentrated  loads,  or  localized 
stresses.   However,  as  noted  and  discussed  in  Appendix  B 
the  general  behavior  of  the  shell  should  be  obtained  first 
and  then  a  follow  up  check  made  for  localized  stresses, 
'The  shell  analysis  should  not  proceed  solely  on  a  series 
of  localized  stress  analyses  throughout  the  structure. 

The  value  of  the  above  equations  is  naturally  dependent 
upon  the  validity  of  the  included  assumptions.   Perhaps  the 
most  questionable  one  of  these  is  the  assumption  of  con- 
stant orthotropic  rigidity  constants  with  variation  in  p„ 
This  simply  states  that  the  structure  is  uniformly  con- 
structed of  the  same  size  structural  elements  around  the 
section.   'This,  of  course,  would  never  be  done  since  the 
loading  varies  with  f.      Figure  26  illustrates  how  the 
actual  coefficients  do  vary  with  $ .      It  can  be  seen  that 


there  are  regions  of  the  shell  where  many  are  constant 
with  0.   For  those  that  vary  slightly,  a  mean  value  might 
be  taken.   The  sharply  varying  ones  are  associated  with 
the  stiffners  in  the  p   direction.   Perhaps  a  reanalysis 
of  the  transverse  strength  requirements  will  reduce 
these  to  more  smoothly  varying  values  with  ft. 

Once  the  equations  are  solved,  the  solutions  will  be 
in  the  form  of  expressions  for  the  deformation  as  a  func- 
tion of  X,  p,  and  the  rigidity  coefficients  at  any  given 
point  I with,  of  course,  the  effect  of  the  injected  boundary 
conditions  being  incorporated  in  the  functional  expression) . 
Once  these  deformations  are  known,  expressions  can  be 
obtained  for  the  strains,  moments,  shears  and  stresses  in 
the  equivalent  orthotropic  shell.   The  rigidity  coefficients 
may  then  be  determined  at  any  point  by  assignment  of  design 
stresses.   These  design  stresses  would,  of  course,  be 
influenced  by  considerations  of  hull  girder  bending  stress 
and  stability  of  the  structure.   After  all  the  rigidity 
coefficients  are  determined  (i.e.,  the  structural  member 
sizes  and  spacings)  they  could  be  fed  back  to  the  original 
expression  for  deformation  along  with  their  variation  in  p. 

A  recycling  of  this  process  until  the  coefficients  do  not 
change  appreciably  will  give  the  final  structural  design. 

A  follow  up  check  for  localized  stresses  and  concentrated 

loadings  will  complete  the  process. 

The  application  of  this  theory  and  procedure  to  a  photo- 
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elastic  scale  model  and  later  a  larger  steel  structural 
model  would  prove  invaluable.   A  program  of  testing, 
modifying  theory  and  reorienting  basic  approaches  would 
most  probably  result  in  new  and  improved  concepts  in 
naval  structural  design. 

A  few  comments  are  in  order  concerning  the  hull 
curvature.   The  reasons  for  selecting  the  analtical 
curvature  (or  radius  of  curvature)  used  in  this  thesis 
are  amply  covered  in  Appendix  A.   However,  for  additional 
research  in  this  area,  NACA  TM  1302,  pp.  8  -  16  is 
strongly  recommended.   The  topic  of  variable  curvature  is 
discussed  and  fourier  polynomials  are  developed  for  the 
radius  of  curvature  expression.   Very  simple  expressions 
result  from  the  analysis  and  they  could  be  used  to  approxi- 
mate, very  closely,  any  given  hull  design.   It  is  quite 
interesting  to  note  that  within  the  scope  of  preliminary 
hull  design,  there  exists  a  design  variable  unused  by  the 
structural  analyst.   This  design  parameter  is  the  miship 
section  coefficient  (the  ratio  of  the  immersed  area  of 

the  midship  section  to  the  product  of  beam,  and  draft,  Ag  )« 

BH 
This  coefficient  defines  the  relative  fullness  at  the 

midship  section  to  a  circumscribing  rectangle.   In  effect 

it  is  a  shape  factor  which  determine  the  shell  section 

profile.   If  the  stress  analyst  had  some  degree  of  control 

over  this  coefficient,  a  shape  or  section  profile  coula  be 

determined  which  would  be  the  most  beneficial  for  the  stress 
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in  the  hull.   The  stress  magnitude  and  distribution  in 
any  shell  structure  is  a  function  of  its  shape  or  curva- 
ture.  The  encouraging  fact  is  that  indeed  the  analyst 
can  play  some  part  in  determining  this  coefficient  to 
the  advantage  of  the  structural  requirements  without 
affecting  the  resistance  of  the  ship.   D.  tf.  Taylor  (8) 
demonstrated  this  fact  in  1908  when  he  conducted  model 
experiments  to  determine  the  influence  of  midship  section 
shape  upon  the  resistance  of  the  ship.   His  tests  were 
definitive  enough  for  him  to  make  the  following  conclusion; 
"for  vessels  of  usual  types  and  of  speeds  in  knots  no 
greater  than  twice  the  square  root  of  length  in  feet,  the 
naval  architect  may  vary  widely  the  midship  section 
fullness  with  no  material  beneficial  or  prejudicial  effect 
upon  speed."   It  must  be  remembered,  of  course,  that  the 
above  statement  applies  to  the  shape  and  coefficient  of 
the  midship  section  for  a  given  area  and  not  to  the  area 
itself.   Taylor  varied  the  midship  section  coefficient 
from  0.7  to  1.1  using  10  models  and  holding  the  midship 
section  area  constant  on  all  tests.   Thus  it  appears  that 
the  structural  analyst  can  accept  a  given  design  and  have 
considerable  freedom  to  determine  his  optlmun  shell  shape 
(for  a  given  midship  area)  within  the  limits  of  not  adverse- 
ly affecting  the    plris>KA/vfc\c  coefficient  or  general 
seaworthiness . 

The  use  of  membrane  theory  for  the  evaluation  of  stress 
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in  the  ships  bottom  section  has  considerable  merit. 
Although  the  stiffened,  ship  structure  will  carry  the 
stresses  according  to  the  distribution  of  its  members," 
membrane  theory  will  predict  what  the  stresses  are  that 
must  be  carried.   There  are  two  basic  limitations  to  the 
theory.   (1)  The  normal  shears,  bending  moments  and 
twisting  moments  are  neglected,  (2)  the  membrane  theory 
becomes  invalid  at  and  adjacent  to  the  boundaries.   These 
limitations  are  not  severe,  however,  since  the  quantities 
of  normal  shear  and  bending  and  twisting  moments  are 
normally  of  a  very  minor  magnitude.   (It  will  be  shown 
shortly  how  this  can  be  determined  quantitatively).   The 
boundaries  can  be  subjected  to  a  separate  discontinuity 
analysis  and  do  not  jeopardize  the  evalution  of  conditions 
over  the  major  portion  of  the  shell.   To  prove  that  the 
membrane  theory  gives  is  a  good  representation  of  the 
situation,  one  can  proceed  in  the  following  manner  (1)  find 
the  displacements  associated  with  the  direct  stresses  pre- 
dicted by  membrane  theory;  (2)   use  these  displacements  in 
the  general  equations  which  incorporate  bending  (equations  10); 
(3)  compute  the  direct  stresses  from  these  bending  theory 
equations  and  compare  them  to  the  membrane  direct  stresses 
to  evaluate  the  accuracy  of  the  membrane  theory. 

Figuresll ,  12,    represent  the  membrane  stress  distri- 
bution in  the  longitudinal  ($=  const)  and  girthwise  (X-  const) 
directions  respectively.   It  must  be  kept  in  mind  that 
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these  membrane  stresses  are  invalid  at  or  adjacent  to 
any  of  the  shell  boundaries  (bulkheads  and  longitudinal 
junctions  with  vertical  ships  side).   A  matching  of 
reaction  forces  at  these  locations  is  necessary  as  well 
as  the  determination  of  secondary  bending  to  satisfy 
compatibility. 

The  most  interesting  feature  of  Figure U  is  the 
presence  of  longitudinal  compression   (-  Nx  stress)  in 
the  shell  midway  between  the  bulkheads.   If  the  ship 
actually  had  a  flat  bottom  the  stress  would  obviously 
be  a  tension.   As  the  bottom  curvature  increases  and 
becomes  a  shell,  the  longitudinal  stress  does  indeed 
become  compression.   The  direct  stress  whose  magnitude 
and  direction  is  the  most  sensitive  to  shell  shape  or 
curvature  is  the  longitudinal  direct  stress,  N  .   This  is 
not  a  conclusion  of  this  thesis  but  rather  a  standard  fact 
of  shell  structure  behavior  (a  fact  which  is   largely 
responsible  for  the  widespread  success  and  popularity 
of  concrete  shell  structures).   This  is  of  interest  to  the 
naval  structural  engineer  since  it  is  this  particular 
stress  which  is  combined  with  the  hull  girder  bending  stress. 
The  recognition  of  this  stress  should  render  ship  structural 
stability  considerations  more  realistic. 

Figure  13  represents  a  plot  of  the  longitudinal  stresses 
as  measured  on  the  destroyer  HMS  ALBUERA  (9)  in  the  hogging 
and  sagging  conditions.   The  test  section  (fr.  62^)  was 
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located,  in  the  midship  region  and  about  halfway  between 
transverse  bulkheads.   Considering  the  observed  stresses 
in  the  hull  above  the  waterline,  a  linear  relationship 
is  readily  established.   When  the  underwater  portion  of 
the  hull  stresses  are  examined  it  can  be  seen  that  they 
exhibit  a  variance  from  the  established  linearity  (the 
linearity  established  the  upper  portion  of  the  hull  must 
be  the  same  as  for  the  lower  since  neither  moment  nor 
section  inertia  has  changed  in  going  to  the  lower  portion) . 
This  departure  is  obviously  due  to  the  effect  of  the  hydro- 
static loading  on  the  underwater  shell  or  hull.   It  will 
be  noted  that  this  effect  augments  the  negative  hogging 
stress  and  reduces  the  positive  sagging  stress.   In  other 
words  the  hydrostatic  loading  on  the  hull  produces  a 
longitudinal  compression  (-Nx).   This  agrees  with  the  theoreti- 
cal plot  of  N  in  Figure  11.   The  longitudinal  compression 
effect  under  hydrostatic  loading  would  be  no  surprise  to  the 
civil  engineer  who  works  with  shells.   It  would,  however,  be 
unexpected  by  a  naval  structures  analyst  who  limits  himself 
to  the  study  of  small  panels  and  predicts  only  a  flat  plate 
tension  effect  due  to  hydrostatic  loading.   While  it  is  quite 
true  that  each  panel  as  such  is  subjected  to  a  tension  pro- 
ducing load,  this  is  only  a  minor  superposition  on  the  more 
general  shell  behavior  which  is  a  compression.   These  local- 
ized stresses  can  be  checked  in  the  later  stages  of  design. 

Figure  12  illustrates  the  distribution  of  direct  stresses 
around  the  girth  (X=  const)  for  two  selected  destroyer  types. 
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The  section  portrayed  is  at  the  left  transverse  bulkhead  of 
the  basic  shell  (X=  -L/2).   This  particular  section  was 
selected  for  purposes  of  obtaining  a  design  shell  thickness 
(and  its  variation  with  0)  since  its  thickness  requirement 
is  greater  than  any  other  intermediate  section  for  any 
given  0   value.   This  was  an  unfortunate  choice,  however,  as 
a  forthcoming  discussion  will  indicate.   The  direct  stresses 
of  Figure  12  were  used  in  calculating  the  equivalent  shell 
thickness  and  hence  extensional  rigidity  constants  (tables  [j.  , 
6,   ).   These  coefficients  are  compared  with  those  actually 
existing  in  the  ships  as  a  result  of  their  construction.   As 
was  previously  stated,  the  membrane  theory  will  give  reason- 
ably good  results  everywhere  except  the  boundaries.   This 
unfortunately  was  the  necessary  region  (trans.  BHD)  to 
analyze  since  it  produced  the  controlling  equivalent  shell 
thickness.   The  evaluation  of  the  direct  stresses  could  have 
been  and  should  have  been  improved  by  using  the  procedure 
indicated  on  page  [j.8  .   Since  the  extensional  rigidity  coeffi- 
cients were  calculated  from  te  which  in  turn  was  calculated 
from  the  direct  stresses;  additional  corrective  action  on 
the  direct  stresses  produced  would  result  in  a  different  and 
better  correlated  set  of  extensional  coefficients.   It  is 
surprising  that  at  this  worst  possible  location  for  membrane 
theory,  the  degree  of  agreement  between  constructed  and  cal- 
culated is  still  reasonably  good.   A  corrective  cycle  on  the 
obtained  membrane  stresses  would  result  in  a  general  increase 
of  varying  degree  for  the  rigidity  coefficients  at  all  values 
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of  0  (since  all  0   values  are  on  the  boundary).   The  differ- 
ence in  relative  variation  of  calculated  from  constructed 
values  in  the  two  ships  selected  is  due  to  several  factors 
which  characterize  each  ship.   It  is  believed,  however, 
that  one  of  the  major  causes  is  the  fact  that  shell  curva- 
ture of  DD  931  (considered  similar  to  DD  952)  resembles  a 
cycloid  to  a  lesser  degree  than  does  DE  1033.   This  sugg- 
ests the  importance  of  selecting  the  most  accurate  analytical 
model  for  hull  curvature.   Other  causes  for  non  similarity 
of  the  ships  plots  are  possible  opposite  departures  from 
the  assumed- conditions  in  the  membrane  equations  (e.g.  ^ 
stress  estimate,  factor  of  safety  used  and  degree  of  fixity 
at  ends).   Figures  15>17  ,   merely  represent  the  comparison 
of  the  calculated  and  constructed  ship  hull  thickness  since 
these  rigidity  coefficients  are  just  a  product  of  a  constant 
and  the  thickness. 

The  approach  used  in  this  thesis  was  not  expected  to 
produce  any  radical  changes  in  structural  design.   It  was 
intended  to  achieve  a  more  rational  design  attitude  by  treat- 
ing a  section  between  bulkheads  as  a  shell  and  then  combining 
these  results  with  the  Overall  ship  bending  stress.   An 
expected  result  was  a  slight  reduction  in  structural  members 
as  determined  by  current  design.   With  due  regard  to  the 
previously  discussed  obstacles  and  handicaps  involved,  the 
plots  of  Figure  11].,  15,16,17,   indicate  an  apparent  trend  in 
this  direction.   The  determination  of  strucutral  members  can- 
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not  be  made  on  extentional  coefficients  alone.   A  recycling 
of  the  membrane  solution  as  noted  on  p   will  give  magnitude 
of  moments  in  the  equivalent  shell.   The  individual  structural 
members  for  the  ship  may  then  be  determined  by  a  combined 
consideration  of  the  extentional  rigidities  required,  the 
moments  in  the  shell  and  the  localized  panel  stress. 

Unfortunately  neither  one  of  the  basic  approaches  used 
in  this  thesis  (i.e.,  shell  bending  theory  by  differential 
equations  and  membrane  theory)  has  been  followed  through 
sufficiently  to  give  definitive  quantitative  conclusions. 
The  investigation,  however,  has  been  quite  profitable  in 
that  the  ships  hull  has" been  analyzed  in  a  manner  more 
appropriate  to  the  geometry  of  the  structure.   Such  an  analy- 
sis can  and,  indeed,  did  reveal  information  unobtainable 
by  an  otherwise  less  realiztic  approach.   Specific  areas  for 
additional  study  and  analysis  have  been  clearly  indicated. 
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V.  CONCLUSIONS 


1.  The  resulting  partial  differential  equations  of  the 
orthotropic  shell  behavior  are  complex  and  even  after 
simplification  are  not  amenable  to  an  analytic  solution. 

2.  A  rational  structural  design  procedure  can  be 
established  by  using  the  partial  differential  equations 
solution  and  the  orthotropic  rigidity  coefficients  (see 
page  4-5" ) . 

3.  The  analytical  expression  for  the  hull  radius  of 
curvature  must  be  uncomplicated  and  yet  representative 
of  the  actual  hull  form  being  analyzed. 

I4..   The  midship  section  coefficient  or  the  midship  sec- 
tion shape  is  a  design  parameter  which  can  be  used  within 
reasonable  limits  to  benefit  the  structural  design  with 
no  prejudicial  effects  on  hull  resistance  or  ship  speed. 
5>.   Subject  to  its  characteristic  limitations,  membrane 
theory  can  be  used  in  the  preliminary  analysis  of  ship 
structures.   Iterative  corrective  procedures  can  be 
used  for  improvement  of  solution  and/or  estimation  of 
existing  error. 

6.   A  longitudinal  compression  results  from  hydrostatic 
loading  on  a  ship  hull  that  possesses  curvature.   The 
magnitude  and  distribution  of  the  compression  around  the 
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girth  of  the  ship  is  sensitive  to  the  hull  shape  or 
variation  of  curvature. 

7.  Useful  extensional  rigidity  coefficients  can  be 
determined  by  membrane  theory  when  iterative  corrections 
are  included. 

8.  Current  design  methods  produce  hulls  which  have  a 
girthwise  bending  stiffness  which  is  five  times  as 
large  as  the  longitudinal  bending  stiffness. 

9.  The  different  and  more  appropriate  approach  used 

in  this  thesis  for  the  analysis  of  ship  structures  with 
hull  curvature  has  produced: 

(a)  A  specific  knowledge  of  items  within  this 
analysis  requiring  further  study  and  refinement. 

(b)  A  more  discerning  attitude  towards  the 
structure  being  designed. 

(c)  Suggested  design  procedures  utilizing  the 
methods  and  techniques  discussed  in  this  study. 
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VI o   REC  OMMENDATI ONS 

1.  Computer  solutions  to  characteristic  partial 
differential  equations  of  the  orthotropic  shell  structure 
should  be  devised,  programmed  and  tested  for  merit  in 
relation  to  current  or  other  proposed  design  concepts. 

2.  Optimization  of  ship  structural  members  and  their 
arrangement  should  be  on  the  basis  of  an  orthotropic 
shell  and  not  on  a  series  of  panels  or  current  design 
approximations . 

3.  A  study  should  be  conducted  to  determine  if  a 
defineable  relationship  exists  between  the  longitudi- 
nal and  transverse  strength  requirements  of  any  given 
design.   The  current  ratio  of  five  to  one  on  girthwise 
flexural  stiffness  to  longitudinal  stiffness  appears 
excessive • 

ij..   Experimental  verification  of  all  developed  theory 
and  techniques  should  proceed  concurrently  to  aid  in 
providing  realistic  results  which  are  of  a  practical 
significance.   Experimentation  could  commence  economically 
by  photoelastic  stress  analyses  of  small  scale  models. 
Tests  on  successively  larger  steel  scale  models  using 
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adequate  instrumentation  should  proceed  if  early  results 
indicate  that  further  testing  can  be  of  value. 

5>.   A  short  study  should  be  made  to  determine  an  analyti- 
cal expression  for  the  radius  of  curvature  of  any  given 
design.   The  expression  must  be  uncomplicated  or  capable 
of  being  simplified  for  the  mathematical  operations 
which  must  be  made  on  it.   (NACA  TM  1302  pp.  8  -  16  is 
recommended  for  this  study) . 

6.  Since  the  mathematical  analysis  of  a  circular 
orthotropic  shell  is  quite  complete,  one  possible  approach 
is  to  break  the  shell  into  a  finite  number  of  segments 
represented  by  constant  radii  of  curvature  (circular), 
make  the  analysis  and  join  the  segments  by  making  them 
compatible  at  all  junctions. 

7.  A  strongly  recommended  topic  for  study  is  the 
variation  of  midship  section  coefficient  or  shape  for  a 
given  midship  section  area  and  the  determination  of  the 
optimum  shape  for  structural  design.   The  only  major 
constraints  will  be  the  maintenance  of  the  designed 
prismatic  coefficient  and  desired  seaworthiness  character- 
istics.  Within  these  constraints,  however,  there  is 
latitude  for  increasing  structural  efficiency. 

8.  A  modified  membrane  theory  should  be  developed  for 
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the  preliminary  analysis  of  ship  structures.   Since  the 
solution  improves  with  iteration,  it  is  quite  amenable 
to  machine  computation. 

9o   A  method  should  be  devised  for  readily  determining 

the  magnitude  of  longitudinal  compressive  stress 
which  is  produced  by  hydrostatic  loading  on  all  hulls 
with  curvature.   This  value  should  then  be  included 
in  any  deliberations  on  structural  stability. 

10,  A  complete  analysis  of  the  stability  of  the 
orthotropic  shell  between  bulkheads  should  be  combined 
with  this  thesis  to  provide  the  complete  ingredients  for 
a  new  concept  of  ship  structural  design. 

11.  A  relieable  solution  to  the  partial  differential 

equation  of  the  orthotropic  shell  will  require  a  know- 
ledge of  the... shell  boundary  conditions.   It  is  apparent 
that  a  study  devoted  soley  to  this  subject  would  be 
of  considerable  value. 
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VII .      APPENDIX 
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A.   HULL  CURVATURE 

The  curvature  of  the  hull  in  the  girthwise  or  circum- 
ferential direction  is  one  of  the  most  significant  factors 
in  this  investigation.   The  author  feels  that  the  magnitude 
and  variation  of  hull  curvature  for  ships  of  the  destroyer 
type  has  a  significant  effect  upon  the  structural  behavior 
and  therefore  should  be  given  full  consideration  in  the 
structural  analysis  of  these  hulls.   The  section  of  the 
ship  considered  was  in  the  midship  section  region.   There 
will  normally  be  no  longitudinal  or  fore  and  aft  hull  curva- 
ture in  this  middle  portion  of  the  ship. 

The  argument  could  be  made  that  the  hull  girthwise 
curvature  is  really  negligible  if  a  basic  panel  (bounded  by 
frames  and  longitudinal  stringers)  is  selected  for  analysis. 
The  weakness  of  this  argument^,  however,  is  that  the  hull's 
general  response  is  not  on  a  panel  basis  but  on  a  much  larger 
and  continuous  region  (i.e.,  the  portion  of  the  hull  structure 
between  transverse  bulkheads)   wherein  any  flat  plate  or 
flat  bottom  assumption  would  be  unrealistic.   Additional 
comments  on  the  relative  size  of  hull  portion  analyzed  will 
be  reserved  until  Appendix  B. 

The  task  of  introducing  the  hull  curvature  of  a  specific 
type  of  ship,  or  of  a  given  design  within  a  type,  is  a  diffi- 
cult one.   The  analytical  expression  of  curvature  selected 
must  be  representative  of  the  design  being  studied.   It  must 
also  be*reasonably  uncomplicated  function  capable  of  being 
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subjected  to  the  necessary  mathematical  operations  without 
excessively  expanding  the  complexity  of  the  resulting 
analysis.   The  representation  of  ships  hulls  by  polynomial 
expressions  is  not  appropriate  since  the  desired  expression 
is  one  giving  the  radius  of  curvature  as  a  function  of  the 
girth  angle  (0)  . 


Figure  18   Cycloid  Profile 
The  analytical  expression  finally  selected  as  satisfying 
the  above  requirements  was  the  equation  of  the  cycloid 
r  =  ij.a  cos  (p      (Figure   )  where  the  depth  of  the  curve 
(draft)  is  2a  and  the  span  (beam)  is  2  IT  a.   ("a"  is  simply 
the  radius  of  the  generating  circle).   This  type  of  curve 
or  section  profile  results  in  a  span  to  depth  ratio  (or  B/H 
ratio)  of  "T  .   The  B/H  ratios  of  the  ship  type  considered 
is  indicated  in  Table  2. 
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DE  1033 

DD  952 

DLG  9 

DLG(N) 

C(x(N) 

CA   134 

B/H 

3.16 

3.13 

3.11+ 

3.02 

3.25 

3. li| 

TABLE  2.   B/H  Ratios  of  Ship  Types  Considered 


The  1/   B/H  ratio  of  the  cycloid  does  not  in  itself 
qualify  this  curve  as  being  a  good  representation  of  the 
destroyer  hull.   The  curvature  and  its  variation  with  jZf 
must  correspond  sufficiently  close  to  the  ship  hull  shape. 
This  is  the  most  significant  factor.   It  must  be  remembered 
that  coincidence  of  the  analytical  expression  and  the  ships 
hull  is  not  of  importance  (i.e.,  the  r,  )6   relationship  should 
not  be  misconstrued  to  be  a  polar  plot).   The  relative 
'equality  of  the  values  of  the  radii  of  curvature  at  given 
0   values  is  the  significant  measure  of  the  expression's 
adequacy. 

The  parametric  equations  of  the  cycloid  in  0   are, 

y  =  a(2jf  4-  Sin2/))  (22a) 

z  =  a(l  +  i22$  (22b) 

The  cycloid  section  can  then  be  fitted  to  any  given  design 
by  letting  a  =  H/2  and  plotting 


y   *    a 


fcos^/z-aUsin 


cos"   /z-a' 


a  / 


obtained  from 


an  e. 


slimination  of  jrf  in  equations  22.   These  plots  have  been 
made  in  figures  20,  21  ,  22  ,  23  ,  along  with  the  plots  of  the 

produced  by  fitting  an  elliptical  section  to  each  actual 


curve 
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hull  shape  selected  for  comparison. 

A  study  of  the  comparative  plots  of  hull  shape  and 
analytical  expressions  reveals  that  the  elliptical  section 
possesses  a  greater  resemblance  to  the  curvature  (or  radius 
of  curvature)  of  the  actual  hull  than  does  the  cycloid  (note 
how  the  DE  1033  curvature  is  virtually  coincident  with  the 
elliptical  section) .   However,  the  analytical  expression 

for  the  radius  of  curvature  of  the  ellipse  is 

2b2 
r  =  f. where  a  and  b  represent  the  major 

a^sin20  4-  b2cos  0 

and  minor  axis  respectively.   This  expression  for  the  radius 
would  not  be  suitable  for  the  mathematical  operations  in- 
volved in  the  shell  theory  of  the  hull.   It  is  also  of 
interest  to  note  how  the  hull  form  of  the  CA  134  is  not 
appropriate  for  this  analysis  since  it  deviates  excessively 
from  curvature  given  by  either  expression.   CG(N)  hull 
form  would  also  be  considered  unsuited  for  the  cycloid  approxi- 
mation. 

The  obvious  disadvantage  of  the  cycloid  section  is  the 
unfortunate  fact  that  the  radius  of  curvature  goes  to  zero 
at  0  -  "/2,      This  is  inherent  in  the  nature  of  the  cycloid 
curve  since  at  the  extreme  points  of  the  span  the  generating 
curve  changes  the  sign  on  the  radius  of  curvature  and,  at 
this  point,  it  passes  through  a  zero  point.   This  is  not  a 
major  disadvantage,  however,  since  it  occurs  at  the  region 
where  the  curved  shell  joins  the  flat  plate  portion  of  the 
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hull  (the  boundary  of  the  shell  section).   This  particular 
region  requires  a  separate  analysis  regardless  of  shell 
theory  employed  (bending  or  membrane  theory) .   In  the 
former  case,  rather  undefineable  boundary  conditions  exist 
along  with  compatibility  obstacles  and  in  the  latter  case 
membrane  theory  never  does  account  for  conditions  at  the 
boundaries.   Therefore,  it  is  of  no  consequence  that  the 
radius  of  curvature  expression  given  by  the  cycloid  is 
invalid  at  the  boundary  of  the  shell  since  the  shell  theories 
are  also  in  question  at  that  region.   A  pertinent  subject 
for  investigation  would  be  the  study  of  the  discontinuity 
stresses  present  at  this  junction*   Results  could  be 
correlated  with  this  thesis  which  purports  to  treat  only 
the  bottom  section  of  the  ship  not  at  or  adjacent  to  the 
boundaries . 
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B.   ORTHOTRQPY  OF  SHIPS  HULL 

As  a  ships  hull  is  fabricated,  the  elements  of 
construction  (shell  plating,  longitudinal  stringers, 
transverse  frames)  are  arranged  in  such  a  manner  that  the 
resulting  structural  elastic  properties  (extensional,  bend- 
ing, shear  and  twisting  rigidities)  are  not  equal  in  all 
directions  as  they  are  for  an  isotropic  structure.   The 
varying  of  elastic  properties  with  direction,  or  the  aniso- 
tropy  of  the  structure,  is  somewhat  simplified  by  the 
symmetry  of  construction  in  the  longitudinal  and  circumferen- 
tial directions.   This  special  case  of  anisotropy  where  the 
elastic  properties  are  different  in  mutually  perpendicular 
directions  is  termed  orthotropy  and  such  structures  are  called 
orthotropic.   It  should  be  noted  at  this  point  that  the  basic 
material  of  construction  is  isotropic.   The  disposition  of 
this  material  creates  the  anisotropy  ana  therefore  the 
structure  can  quite  correctly  be  called  constructionally 

orthotropic. 

The  reaction  of  a  structure  to  any  given  loading  will 
be  a  function  of  its  elastic  properties.   It  would  seem 
tjien  that  a  realistic  analysis  could  be  made  if  the  structures 
elastic  properties  and  their  change  with  direction  are  recog- 
nized.  This  can  be  accomplished  by  considering  an  "equiva- 
lent" structure  without  stringers  and  frames  but  whose  elastic 
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properties  and  general  behavior  under  load  resembles 
sufficiently  the  original  structure. 

The  advantage  of  this  approach  is  that  the  continuity 
of  the  structures  behavior  is  retained  and  the  resulting 
analysis  is  superior  to  one  in  which  conservative  attempts 
are  made  to  account  for  the  complex  interactions  existing 
between  structural  elements  such  as  frames,  longitudinals 
ahd  shell  plating.   A  critical  question  arises  at  this 
point  -  what  is  the  characteristic  mode  of  deformation  of 
tjie  stiffened  hull?  Does  it  deform  under  load  to  form  a 
smooth  V^aH  wave  length  between  bulkheads  as  in  Figure  2lj.a 
or  is  all  the  deformation  confined  to  the  shell  surface 
between  stiffening  members  as  in  Figure  2l4_b  ? 


^FT1  r  A"J 
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14cl. 


t  I 
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Figo  2k     Stiffened  Hull  Modes  of  Deformation 
The  answer  is  clearly  a  function  of  the  number  and  size  of 

the  transverse  frames.   Lundgren  (7)  is  satisfied  to  accept 
the  former  version  of  deformation  (Figure  2lj.a) ,  whenever  the 
number  of  transverse  frames  between  bulkheads  equals  or 
exceeds  four.   There  are  panel  advocates,  however,  who 
insist  on  the  latter  mode  (Figure  21+b)  as  being  the  predomi- 
nant one  which  accounts  for  most  of  the  strain  energy  of  de- 
formation regardless  of  the  number  of  frames. 

For  the  destroyer  type  having  approximately  i\.   frames 
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between  bulkheads  (in  the  middle  half  of  the  ship),  the 
author  feels  that  the  normal  mode  of  deformation  is  probably 
a  combination  of  smooth  half  wave  length  deformation  with 
smaller  panel  deformation  superposed  on  it  as  in  Figure  25  . 


^4b 
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Fig.  25  Combined  Modes  of  Stiffened  Hull  Deformation 
Perhaps  many  of  those  who  can  accept  only  panel  deformation 

are  misled  by  the  physical  appearance  of  ships  hulls  which 
have  the  so  called  "hungry  horse"  look  or  a  permanent  deforma- 
tion on  all  panels  making  the  shell  appear  as  though  it  were 
shrunk  fit  onto  the  framing  system.   These  effects  could 
certainly  have  originated  from  shrinking  of  the  shell  around 
faying  surfaces  welded  to  the  hull  or  possible  a  dynamic 
loading  effect  of  the  ocean  waves. 

The  deformation  of  Figure2£  (which  is  the  authors 
version  of  the  situation)  illustrates  that  continuity  of  the 
structure  between  bulkheads  must  be  considered  in  analysis. 
The  extent  of  this  significant  structural  section  also  requires 
a  full  consideration  of  shell  curvature.   One  handicap  of 
the  gross  section  approach  which  is  also  obvious  from  Figure 
is  that  the  localized  stresses  created  by  the  superposed 
panel  deformation  would  not  be  immediately  determined.   A 
final  check  of  the  solution  should  incorporate  a  method  of 
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accounting  for  these  localized  stresses. 

The  orthotropic  rigidity  coefficients  (equation  II4.)  have 
been  obtained  and  defined  by  the  superposition  of  the  elastic 
laws  of  the  isotropic  shell  and  the  mating  gridwork  frame. 
A  better  appreciation  for  these  may  be  achieved  by  an  examina- 
tion of  these  coefficients  on  an  actual  hull.  Figure  26 
illustrates  the  magnitude  and  variation  of  the  coefficients 
for  the  DD  931  class  destroyer.   The  data  for  these  plots 
is  in  Appendix  D   . 
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Fig.  26   Orthotroplc    \i 
constructed  on  DD 


Rigidity  Coefficient^  as 
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C.   BOUNDARY  CONDITIONS  FOR  USE  wilTH 
PARTIAL  DIFFERENTIAL  EQUATIONS  OP  SHELL 

The  solution  of  the  developed  partial  differential 
equations  {IS  )    will  eventually  require  the  stipulation 
of  appropriate  boundary  conditions.   These  conditions  will 
consist  of  a  rational  approximation  of  (a)  the  degrees  of 
fixity  existing  at  the  bulkheads  and  at  the  junction  of 
the  curved  shell  and  flat  ships  side  and,  (b)  the  magnitude 
and  distribution  of  the  loading  which  exists  at  the  afore- 
mentioned boundaries. 

Fixity; 

The  degree  of  fixity  which  exists  at  the  boundaries 

of  the  shell  could  indeed  be  the  subject  for  a  separate 
and  absorbing  analysis.   Vedeler  (1)  suggests  a  method 
for  accomplishing  this.   It  would  be  of  questionable  value 
however,  to  pursue  this  subject  in  any  depth  in  this  thesis 
since  the  immediate  problem  of  de terming  a  method  of 
solution  to  the  differential  equation  is  the  major  obstacle. 
For  a  first  approximation  then,,  the  following  conditions 
are  postulated.   At  the  junction  of  the  shell  and  the  bulk- 
head, the  fixity  can  be  considered  to  be  100$  clamped.   The 
presumtion  is  a  natural  result  when  one  considers  that  the 
structure  is  a  shell  continuous  over  supports  (bulkheads) 
with  continuous  and  equal  loading  on  each  side  of  these 
supports.   An  examination  of  the  bulkhead  Nx  stress 
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distribution  in  Figure  11  would  also  lead  one  to  assume 
that  a  full  clamping  condition  prevails.   The  fixity 
condition  along  the  longitudinal  boundary  is  a  bit  more 
elusive.   Once  more,  as  a  first  approximation,  it  will  be 
assumed  as  being  simply  supported  and  capable  of  movement 
in  the  athwartship  direction.   The  condition  can  be 
represented  by  Figure  27  • 


Figure  27  Degrees  of  Fixity  at  Shell  Boundaries 
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Edge  Loading; 

To  begin  the  estimate  of  edge  loading  and  its  dis- 
tribution, consider  the  bottom  section  and  the  total 

ON  '* 

vertical  force^due  to  the  water  (q)  and  the  internal 
loading  (w) . 

q  =  q0  4  Y  z 

where  qQ  -  y hQ   ,  hQ  =  a  design 

head  for  protection  of  freeboard 
plating. 

w  '  const,  across  section  (^/n1-) 

qv  -  (%+  X   z)cos0 

qv  =  V"(h0-\-  z)cos0 

r  *  2H  cos0 

dS  =-  rLd0 

z  =  H/2  (1  4-  cos2J()) 


Fig.  28  Shell  Loading 

total  vertical  force  on  hull  (\J.  \ 

t   =      /Uv  -  w)dS 
=  2     J(qv  -  w)rLd|Z> 


V 


/% 


=  2     J  (    y  hQ   cos0   -    Yz   cos$  -  w)2H  cos0  Ld0f 
-  k-HL^  (hQ    cos2jZ$  -|-  z    cos20   -  w   cos0)    d0 

=  ij.HL 


* 


Vt   -     HL 


Pcos20    +H/2(cos20  4-  cos20cos20j    -  w  cos^l     d(j) 
r       ™ 

if    ffhnf   3f  H-i[Wj 
4  ft 
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An  alternate  expression  for  V.  could  have  been, 


vt  - 


(q  -  wn  -  wt)  dSr 


where  wn,  w^  are  normal  and  tangential  components  of  the 
internal  loading  w.   dSzis  an  increment  of  surface  area 
projection  on  the  horizontal  plane. 
Substitution  of  wn  ■=  w  cosO 

wt  -s.  w  sinO 

dSz  —  f   L  cos0  d0 
and  integration  will  give  same  results  obtained  previously 
for  V£. 

Having  now  obtained  the  total  upward  force,  a 
determination  must  now  be  made  of  the  distribution  of  the 
reaction  forces  which  support  V-^  and  act  on  the  shell 
boundaries.   This  is  a  statically  indeterminate  problem. 
In  order  to  circumvent  this  obstacle  and  avoid  an  auxiliary 
analysis  ^ whose  accuracy  would  be  greater  than  justified^, 
the  following  approach  is  taken. 


i 
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Fig.  29  Analogy  of  Flat  Plate  Prismatic 
Loading  to  Hull  Loading- 

-78- 


Timoshenko  (10  describes  the  behavior  of  a  simply 
supported  rectangular  flat  plate  under  a  triangular 
prism  distribution  of  loading,  Figure  29.   The  portion 
of  this  reference  which  is  of  interest  is  the  nature  of 
the  reactive  forces  along  the  simply  supported  edges. 
Indeed  it  could  be  said  that  the  vertical  component  of 
loading  on  the  ships  hull  is  a  triangular  prism  distri- 
bution (with  max  value  at  the  keel  J  if  the  hull  were 
developed  out  into  a  flate  plate,  Figure  29.   Because 
the  ships  loading  is  not  on  a  flat  plate  but  on  a  cylin- 
drical hull,  no  further  resemblance  or  analogy  will  be 
made  other  than  a  consideration  of  the  relative  distribu- 
tion of  the  reactive  forces  and  their  relationship  to  the 
aspect  ratio  of  the  simple  supports.   Timoshenko  tabulates 
VY  max  '  Vx  max  for  various  aspect  ratio  of  simple  support 
boundaries.   For  the  ship  types  being  considered  in  this 
analysis,  the  pertinent  ratio  of  B/L  is  approximately  1.3 • 
This  value  gives 

vy  max     2.73  x  aspect  ratio 
^x  max 

or 


2.73  L/B 


VB 
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Therefore, 


ft 


2LVL  +2BVB     =      VT 


2LVL  +    2   BVB    =    HL  V    (iThQ  i-    3A^    H  -   l±w  j 


and, 

VL  2.7    L/B 

Solving  for  the  reaction  forces  gives, 


rL  =  1.35  BHL2  Y  (  ir  h0  4  3JT  H  -  i|.  w J    |t. 
2.7  L2  4-B2  l+         V 


V   =  Q.5  BHL  Y  (   <  h0  +  2J[  H  -  I4.WJ 


2.7  L2+B2 


Pig.  30  Edge  Loading  on  Shell  Boundaries 
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Rl     ftlGI^ITY  COEFFICIENTS  DATA 

The  orthotropic  rigidity  coefficients  of  equation  (IJ4.) 
can  be  written  in  the  following  form  for    S)-  0.3. 
Extensional  coefficients 

D_    Lit    A 

JL   -  V     — 

E  b1 


DJ.  ■■ 

l.lt  ^  A^ 

E 

b2    > 

2^ 

0.33t 

E 

^ 

0.38t 

E 

Flexural 

coefficients 

K 

X 

O091^t3    Ix 

2 
Vx 

E 

.          |.  —  — 

b2 

2.  -■ 

E 

4-  - 

bi 

S.i    . 

.06Ip.t3 

4 

.38  jx 

E  b2 
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K0x  ._  .061^-lt3  x     .38  J0 

>1 


E  b- 


Kv  z    .0275t3 

E 

Rigidity  moments 
Sx    Axcx 


E       b2 


E     bn 
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TABLE  3   DP  931   RIGIDITY  COEFFICIENTS  DATA 


0 

0° 

10° 

30° 

1^5° 

72° 

90° 

(in)    t 

.75 

.75 

.625 

.i+375 

.3125 

.375 

(in*)   Ax 

19.8 

2.53 

2.53 

2.53 

2.014- 

2.QI4. 

(in)   b2 

38.it- 

38.4 

38.4 

38.14. 

28.8 

28.8 

(in)    Cx 

17.03 

6.73 

6.66 

6.57 

k.3k 

14-37 

(H*)    Ix 

IIOJ4.. 

26.32 

26.32 

26.32 

7.61 

7.61 

(in^-)    Jx 

12x5 

27.37 

27.37 

27.37 

8.81 

8.81 

(in2)   A0 

13.5 

12.75 

8.7 

8.7 

5.61 

5.21+ 

(in)    C0 

17.83 

16.38 

12.66 

12.57 

9.66 

9.1*0 

(in1*-)   10 

859.5 

650.5 

328.5 

328.5 

III4..9 

107.8 

(in^)   J0 

931.5 

722.5 

3¥i«5 

3li4.5 

121.7 

111.8 

(in)        b     =     8I4.  =■      const  0 
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TABLE  k     DP  931   CALCULATED  EXTEflSlOML  HiGIDITY 
COEFFICIENTS  AT  BULKHEAD 


i 

0 

0 

0 

10 

0 
30 

k$' 

0 
60 

0 
72 

80° 

Nx 

+13936 

+13341 

+  9150 

+  3200 

-3730 

+  I4.80O 

1^8000 

V 

•35805 

-314-360 

-214.950 

-15,380* 

-7360 

-  B£2^ 

"  m? 

u 

0 

10753 

27910 

33920 

33800 

32690 

38850 

ri 

16600 

16220 

13450 

10^70 

7300 

5310 

14480 

*. 

1.215 

1.08 

0.58 

0.36 

0.308 

0.293 

0.88 

E  "  E 

1.339 

1.19 

.638 

.396 

.339 

.322 

.969 

E 

o.k 

0.356 

0.191 

0.119 

0.102 

.096 

0.29 

x0 
E 

0.461 

0.I4.I 

0.22 

0.137 

0.117 

0.111 

0.331+ 
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TABLE  5   DP  931  RIGIDITY  COEFFICIENTS  (AcIuclL) 


t 

0° 

10° 

30° 

hS° 

72° 

90° 

s 

1.3W 

0.891 

0.7514- 

0.514-7 

O.I4J.5 

O.I4-8I4. 

D0 

0.986 

0.977 

0.792 

0.585 

0.1*11 

0.14.75 

E 

O.2I4.8 

O.2I4.8 

0.206 

o.ilj4 

0.103 

O.I2I4. 

£ 

0.285 

0.285 

0.238 

0.166 

0.119 

0.11|.3 

E 

178.8 

3.70 

3.614- 

3.5I4- 

1.60 

1.62 

K0 

E 

61.  k$ 

UB.39 

20.63 

20.36 

7.56 

6.7 

2 

0.0115 

0.0115 

0.0067 

0.0023 

0.0008 

0.0015 

E 

12.03 

0.298 

0.287 

0.276 

0.118 

0.120 

E 

1+.237 

3.297 

1.576 

1.565 

0.552 

0.539 

Sx 

E 

8.8 

0.1* 

0.144 

O.I4.3I4- 

0.308 

0.310 

E 

2.87 

2.14-8 

1.32 

1.31 

O.6I4- 

0.58 
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TABLE  6  DE  1033  CALCULATED  EXTENSIONAL  RIGIDITY 
COEFFICIENTS  AT  BULKHEAD 


0 

0° 

10° 

30° 

45°  ■ 

60* 

0 

72 

80° 

db    H 
it)    x 

4  7725 

.7414 

t5201 

+2180 

-750 

+  8000 

+110,000 

N* 

-2^,580 

-23650 

-17330 

-10880 

-5375 

-2463 

-1162 

0 

6463 

16910 

20,770 

21,250 

21,600 

27,200 

(-hi)    fi 

15000 

14680 

12150 

9450 

6600 

4800 

4050 

M*e 

0.538 

0.495 

0.321 

0.240 

0.198 

0.202 

0.635 

£        £ 

0.592 

0.545 

0.354 

0.264 

0.218 

0.222 

0.7 

£ 

0.1775 

0.103 

0.106 

0.0792 

0.0654 

0.666 

.21 

^X0 

E 

0.204 

oa88 

0.122 

0.091 

0.0753 

.0768 

0.242 
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TABLE  7   DE  1033   RIGIDITY  COEFFICIENTS  (Ac^o-l) 


0 

0 

0 

10 

30 

o 

45 

0 

60 

0 

72 

0 

80 

t 

0.5 

0.5 

0.4375 

0.4375 

0.375 

0.375 

0.3125 

** 

24.71 

2.81 

2.49 

2.20 

2.20 

2.20 

2.20 

b2 

26.9 

26.3 

27.5 

25.6 

26.3 

28.1 

26.2 

v 

8.5 

7.3 

4.63 

4.57 

1|.27 

4-03 

3.91 

E 

1.47 

0.657 

0.571 

0.566 

0.497 

O.489 

0.428 

St 

E 

0.644 

0.631 

0.531 

0.53 

0.46 

0.458 

0.387 

E 

0.165 

0.165 

0.144 

0.144 

0d24 

0.124 

0.103 

Dx<3 

E 

0.19 

0.19 

0.166 

0.166 

0.142 

0.142 

0.119 

b1  -      90    = 

-     const. 

i 

| 

1 
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